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SATELLITE BROADCASTING SERVICE AREAS 



SUMMARY 

The coverage area provided by a transmitter on a geostationary satellite 
emitting a right-circular conical beam is mainly considered. For sufficiently 
small beamwidth, the service area is nearly an ellipse in the tangent plane to the 
Earth at the target point, and the axial lengths and directions are calculated. In 
general, the service-area boundary is best specified by plotting latitude against 
longitude, and formulae for specifying this boundary for arbitrary beamwidth are 
given and discussed. In particular, the 'tulip' shape of some of the diagrams of 
Fig. 2 (reproduced from reference 1) is explained. 



1. INTRODUCTION 

If a satellite is to be geostationary, it must be 
above a point on the Equator, and at a distance of 
about 42,200 km from the centre of the Earth. The 
beam emanating from the satellite will generally be 
assumed to be a right-circular cone (the case of a 
general conical beam isdiscussed in outline inSection 
5). This cone will meet the Earth's surface in a 
curve, as illustrated in Fig. 1. For sufficiently small 
beamwidth, the curve is approximately an ellipse in 
the tangent plane to the Earth at the target point. For 
larger beamwidth, the curve will cease to lie in one 
plane, although it will remain a closed curve unless 
the beamwidth is so great that part of the cone of 
radiation from the satellite lies outside the tangent 
cone from the satellite to the surface of the Earth. 




Fig. 1 ~ Illustration of a typical service area (for 
narrow beamwidthl 



In general, it is more convenient to specify the 
boundary of the service area by plotting latitude 
against longitude. Fig. 2 {reproduced from reference 
1) indicates the general nature of such curves. 
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Fig. 2 ~ Coverage of Europe with the satellite at 8°E, Paris as the target point and varying beamwidth 
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The object of this report is toprovide formuiae by 
means of which the boundary of the service area can 
be determined whatever the beamwidth. The notation, 
based mainly upon the geometry of Fig. 3, Is explain- 
ed in Section 2. The case when the beamwidth is 
smail is considered In Section 3; axial ratios and 
directions are tabulated for ten-degree intervals of the 
latitude of the target point, and of the difference 
between the longitudes of the target point and the 
subsatellite point. 

If latitude is plotted against longitude, it can be 
shown that the small-beamwidth service-area boun- 
daries remain ellipses, but as the beamwidth increases, 
they become non-elliptical closed curves until the 
cone of radiation from the satellite meets the tangent 
cone from the satellite to the Earth. For still greater 
beamwidth, we have the 'tulip' shape of some of the 
curves In Fig. 2. To determine a typical point of such 
a curve, it Is necessary to find first the coordinates of 
a geometrically-specified point 0^, on the boundary 
of the service area and in the plane through the target 
point Pi normal to the line SPi joining Pi to the 
satellite S. Thereafter we have to determine whether 
the line SQ„ meets the Earth's surface, and, If so, 
the coordinates of P^i the point of intersection. 
Formulae determining the coordinates of 0^. and the 
coordinates of P^, in terms of those of Q^, are given 
in Section 4. Section 6 gives conclusions. 



2. NOTATION AND THE GEOMETRICAL CONSIDER- 
ATIONS UNDERLYING IT 

Most of the points which are relevant to the deter- 
mination of the service area are in the plane SCPi 
(the plane of the paper in Fig. 3), where S is the 
satellite, Pi the target point and C the centre of the 
Earth, The coordinate axes (not shown in Fig. 3) are 
such that Cz is the axis of the Earth's rotation, and 
Cx is the intersection of the Equator and the Green- 
wich meridian; for a right-handed system of axes, lati- 
tude North, and longitude East of Greenwich, must be 
taken as positive. We shall use the following sym- 
bols:- 

C The centre of the Earth (and origin of co- 

ordinates) 
D The distance CS (taken as 26,300 miles 

(42,200km)'in numerical work) 
Ki The point on the Equator due South of Pi 

Li.Mi Arbitrary points In the plane SCPi, such that 
LiPiMi is the tangent to the Earth at Pi in 
that plane. 
Pi The target point 

P„ {n i= 1)A point where the line SQ^, meets the Earth's 

surface 
Q2, Q.i Points in the plane SCPi, in the plane 
through Pi normal to SPi and on the boun- 
dary of the service area 
Q„ int- 1,2,3) 

Any other point on the boundary of the ser- 
vice area and inthe plane through Pi normal 
to SPi 




The distance Q^Pi (= SPitanVia) 

The radius of the Earth (taken as 3960 miles 

(6375km) in numerical work) 

The subsatellite point 

The position of the satellite {co-ordinates 

(Dcosc^Ss, Dsin05, 0)) 

(= Rcos^„cos0„) X — coordinate of P„ 

(= Rcos^^sin^^) y — coordinate of P^ 

(= Rsln0„) z — coordinate of P„ 

Coordinates of Q„ 

The beamwidth 

The angle between CPi produced and PiQ^ 

The latitude of F= 

The angle Q^PiQa 

The longitude of P„ 

The longitude of s 

The angle SCPi 

The elevation of S as seen from Pi 

The bearing (West or East of South) of s from 

Pi 




Fig. 3 - Relevant points associated with the plane 

formed by the satellite, the target point and the centre 

of the Earth 



3. AXIAL LENGTHS AND DIRECTIONS FOR THE 
NARROW-BEAM CASE 

From the geometry of Fig. 3 the following rela- 
tions can be established:— 

Dcosi/zi-R 

siny^=^ ; cosy^ = Dsint/'i/SPi (1) 

SPi 

cosy^ - cosy^ cos^ = Dsini^i cos^/SPi (2) 



SPi =D^ +R''-2DRcosiP, 

COSl/'i =0056*1 COS(0i~(^s 



(3) 
(4) 



and from the spherical triangle PisKi, right-angled 
at Ki 



sinciii =r^ sin((i.-cii)/sin(/'i 



(5) 



For the narrow-beam case, we can regard the 
beam from S as a right-circular cylinder of radius r, 
cut by the tangent plane to the Earth (in which the 
elliptical service-area boundary lies) at an angle 
(n-/2 - yj with the plane through Pi normal toSPi. 
Hence the major axis of the ellipse is in the direction 
LiPiMi and has length rcosecy^. and the minor axis 



(in the direction through Pi perpendicular to the 
plane SCPi) has length r. The ratio cosecy^ of the 
axial lengths can be deduced directly from the first 
of (1), and is tabulated in Table 1, while the direction 
of the major axis is determined by coi in Equation (5) 
and is tabulated in Table 2. 



TABLE 1 



Axial Ratio 



^\ ^1 


0'^ 


10° 


20° 


30° 


40° 


50° 


60° 


70° 


80° 


'^s-'KX. 




















0° 


1 


1.021 


1-090 


1-220 


1-446 


1.850 


2-673 


5-009 


42-68 


10° 


1.021 


1.043 


1-113 


1-246 


1-477 


1-891 


2-736 


5-153 


48-21 


20° 


1.090 


1.113 


1-188 


1-330 


1-577 


2-022 


2-940 


5-633 


78-25 


30° 


1-220 


1.246 


1-330 


1-489 


1-769 


2-277 


3-344 


6-634 


. - 


40° 


1.446 


1.477 


1-577 


1-769 


2-108 


2-736 


4-098 


8-718 


- 


50° 


1.850 


1.891 


2-022 


2-277 


2-736 


3-609 


5.633 


14-12 


- 


60° 


2-673 


2-736 


2-940 


3-344 


4-098 


5-633 


9-789 


48-21 


- 


70° 


5-009 


5-153 


5-633 


6-634 


8-718 


14-12 


48-21 


- 


- 


80° 


42-68 


48-21 


78-25 


- 


- 


- 


- 


- 


- 



TABLE 2 

Direction of Major Axis 
Angle (oii) from aNorthSouth Line is Tabulated 



^\ ^1 


0° 


10° 


20° 


30° 


40° 


50° 


60° 


70° 


80° 


^s-^lX^ 




















0° 


- 


0° 


0° 


0° 


0° 


0° 


0° 


0° 


0° 


10° 


90° 


45°27' 


27°16V2' 


19°25V2 ' 


15°20' 


12°58' 


11^0V2' 


10^8' 


10°9 ' 


20° 


90° 


64° 30 ' 


46°47 ' 


36°3 ' 


29°3r 


25° 25 ' 


22°48 ' 


21°10' 


20°17' 


30° 


90° 


73°16' 


59° 22 ' 


49°6 ' 


41°56' 


37°0' 


33°41V2 ' 


31°34' 


- 


40° 


90° 


78° 181/2' 


67° 49 ' 


59°13' 


52°33 ' 


47° 36 ' 


44°6' 


41°46' 


- 


50° 


90° 


81°42V3' 


73° 59 ' 


67°14V2' 


61°39V3' 


57° 16' 


54°0' 


51°45' 


- 


60° 


90° 


84° 16' 


78°50 ' 


73°54' 


69°38V2 ' 


66°9' 


63° 26 ' 


61°31 ' 


- 


70° 


90° 


86°22V2 ' 


82° 54 ' 


79°41 ' 


76°50 ' 


74° 25 ' 


72°30 ' 


- 


- 


80° 


90° 


88°14' 


86°33 ' 


- 


- 


- 


- 


- 






4. DETERMINATION OF A TYPICAL POINTOF THE 
BOUNDARY OF THE SERVICE AREA IN THE 
GENERAL CASE 

When the beamwidth is not necessarily small, we 
seek to determine the latitude 0^ and longitude 0^, of 
a typical point P^ of the boundary of the service area, 
and for this purpose it is convenient first to find the 
coordinates of the point O^ in which SP„ meets the 
plane through Pi normal to SPj,. We assume that 0^ 
is specified in terms of r (which is determined by the 
beamwidth and the position of the target point Pi, and 
is thus a known constant) and the angle ^^ or O^PjOs. 
Then since PiQ^j = r 

2 2 2 2 

X, ■,Y,+Z,=r (6) 

Since PiQ^ is perpendicular to SPi, 
(Dcos05-Rcos{SiCOSQbi)X^ + 

+ (Dsin(/>s-Rcos^isln(/)i)y„-Rsin(9iZ^ =0 (7) 

and since the angle between CPi produced and PiQ,] 
is y„ specified by equation (2), it can also be shown 
that 

cosi9iCOS(^iX^ + cosf^iSinc^iVn + sinOiZ^ 

= rDsiniAiCos4/SP, (8) 

In general, (that is, if (f>i / '^s) (^)> ('') ^"^ (^' ^^^ ^^ 
reduced to 

2 2 r" T 

Z„ sin ^1 -2rZ,jCoSif„sin^isini/-'i[_D-Rcosi/'iJ/SPi 

?. i 2 2 2 ^ 

+ r <sin i9i-sin ly/isin ^^ 
- (R^sin^f?iSinViCos^4^SPi^) > = 
and X^, Y„ are then obtained from 



(9) 



;10) 



where 



X^ sin(0i - 0^) ^ B sint^^i - A sin^^ 
y„ stn(0i - q!»s) - A coscf)^ - B costf^i j 

A = j(rDsiniAiCos4,)/SPi^Z^sin^][/cosf?i} 

B= (rKsiniAiCosf^)/SPi ( 

(11) 

Given r, i^ [as well as 9-^, (0i--<jSg) and hence i/fil 

(9) gives two values of Z„ as a multiple of r, and then 

(10) and (11) give the corresponding values of X^ and 
y^ as multiples of r. Equations (9), (10) and (11) are 
unaltered if ^„ is replaced by -^„; this explains why 
a quadratic equation for Z^ is to be expected. If 
01 = 0s, the corresponding solution is explicit, name- 
ly 

X„ = r[+ sin0isin^„ + (K/SPi) sin(9i oos0i cos^^J 

^n = '■[+cos0isin^,i 1 (/?/SPi) sini9i sin0i cosifj 

Z^ -:: rcos4 {D - ficos^i)/Sp3 

(12) 



At this stage we are in a position to determine 
explicitly the coordinates of 0,^, which are 

(xi + x^, >'i + y^, zi + z^), 

when r and ^„ are given; it remains to determine the 
latitude and longitude of the point P„ in which SQ^ 
meets the Earth's surface. It is sufficient to deter- 
mine the coordinates (%^, y,j, z^) of P^; it can be 
shown that these satisfy the equations 



x^-%1 -X^ 



yn-yi-yn 



Dcos0s-xi-X„ Dsin0s-yi-Y^ 

Z 1 + ^ri ~'^r 



2i +z, 



= ^y. say 



:13) 



where 



T]'' [spf +r^ - (2rKDsini/'iCos4,/SPi)] 



-~277[DRcosi/fi +(rRDsin(/'aCos4,/SPi)-R -r J + r =0 

(14) 

Equation (14) (of which only the smaller root is 
relevant) only involves ^^, r and known quantities, 
so T] can be determined without knowing the position 
of Q^ explicitly. Once q is known, the last of equa- 
tions (13) gives z^ and hence ^^ directly. Either of 
the remaining equations (13) then gives x^ or y^ and 
hence 0^!- 

We are therefore now able to plot the latitude 
6^ against the longitude 0„ in the general case, using 
a linear scale for both quantities. The geographical 
significance of the results is best appreciated by ap- 
propriately distorting the map of the world in the 
relevant region, since It is only major geographical 
features with which we are concerned. 

If the beamwidth a (and therefore r) is sufficiently 
small, the line SO2 in Fig. 3 will meet the Earth's 
surface at a point P2 (not shown in Fig. 3) on the 
great circle sPi produced, but if the beamwidth is too 
large, the line SQ2 will fail to meet the Earth. The 
critical situation occurs when the line SQ2P2 tpuches 
the Earth at P2, so that the angle P2SC is sin" [R/D) 
or 8^40'. If ang'e P2SC exceeds this value, there 
must be a value of ^^ such that angle P^SC is 

sin'^(E/D) 

instead, and the point P^ represents the theoretical 
extremity of the 'tulip' in such a case. At this point 
P„, the elevation of the satellite is zero, so that 
service, though theoretically possible, will be poor. 
In practice, the service area should be regarded as 
terminated when the elevation of the satellite has a 
minimum value of at least a few degrees: the 'tulip*- 
shaped curves of Fig. 2 are for this reason terminated 
before the theoretical extremity is reached. 



5. MODIFICATION FOR A CONICAL BEAM OF EL- 
LIPTICAL CROSS-SECTION 

If the beam emitted from the satellite is an arb- 
itrary cone, meeting the plane through Pi normal to 
SPi in an ellipse with centre Pi, the only difference 
will be that the distance r to a point on the service- 
area boundary in the plane through Pi normal to SPi 
will no longer be constant and equal to SPitanVaa, 
but will instead be given by a formula of the form 

1 1 2 1 2 

-2=-2C0S {4-^o)+-2Sin (4~^o) (15) 
r ri r^ 

where r-i = SPitanViai Is the distance to a point of the 
service-area boundary in the direction of maximum 
beamwidth Viai, ^o is the value of ^n '" ^'^'s direc- 
tion, and r2 = SPitanV2a2 is the distance to a point of 
the service-area boundary in the perpendicular direc- 
tion of minimum beamwidth Via^. In equation (15) there 
are three distinct essential parameters, r^, r^ and ^o. 
whereas when ri = T'z the value of ^o is immaterial 
and we only have one parameter. 

Qualitively, we can expect the same general 
features as before; nearly elliptical service areas 
for sufficiently small beamwidth, and 'tulip'-shaped 
service areas for larger beamwidth. But as the orienta- 
tion of the beam now has to be considered, there Is no 
longer necessarily symmetry of the service area boun- 
daries with respect to the plane SPiC of Fig. 1, and, 
as already noted, there are three significant para- 
meters to consider instead of only one. No detailed 
quantitative investigation will therefore be under- 
taken here, although ail the mathematical machinery 
for such an investigation is available in the formulae 
already derived if a practical requirement should arise. 



6. CONCLUSIONS 

For a right-circular beam of sufficiently small 
beamwidth, the service-area boundary is approximate- 
ly an ellipse In the tangent plane to the Earth, centred 
at the point towards which the beam is directed. The 
axial lengths and directions of this ellipse are deter- 
mined by simple formulae based on the geometry of 
Fig. 3. For larger beamwidths the service area boun- 
daries have the 'tulip* shape so noticeable in Fig. 2 
(taken from Reference 1); formulae are given for deter- 
mining the position of an arbitrary point of such a 
service area, and in particular, for determining the 
theoretical extremity of the 'tulip'. The case of a 
conical beam which is not right-circular is qualitative- 
ly similar and could be fully investigated by means of 
the formulae given here, but has not been attempted 
because it is qualitatively much more complicated, 
since three significant parameters are involved in- 
stead of one. 



7. REFERENCE 

1. BENoIt, a., GODFROID, H. and KNYPERS, P. 
1968. Distribution of television by satellite, with 
special reference to the size of the Earth-station 
aerials. E.B.U. Rev., 1968, Part A, No. 110, 
pp. 162-172. 



JMP 



Printed by BBC Research Department, Kingswood Warren, Tadworlh, Surrey 



read 



where 



CORRIGENDA 

RESEARCH DEPARTMENT ~ BRITISH BROADCASTING CORPORATION 

Research Department Report No. 1970/3 
and 
BBC Engineering, April 1970, No. 82, p. 28 

SATELLITE BROADCASTING SERVICE AREAS 

The author regrets that he has recently found out that Equation (14) of the original text is wrong, and should 

7?=(SP^+r2)+T7[(D^ -SP^)-(/-2 +Ca„^)] +CQ„2 -R'' =0 



CQ„2 =i^2 +r=' + {2rRDsm \p ^ cos ^„)/SP^ 

It may be useful to note that Equation (9) also has an explicit solution 
Zfj=rs\n6^ cos^^[D- R cos i^/^] /(SPj sin \p^)±r cos 6^ sin ^„ sin(0^ ~0g}/sin \p ^ 
and this can be used to simplify the formulae (10) forX^ and Y^ in a form that covers ail cases. 
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